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We explore the phase structure of dense magnetized quark matter when a repulsive vector inter¬ 
action, parametrized by Gv^ is present. Our results show that for a given magnetic field intensity 
(B) one may find a value oi Gy for which quark matter may coexist at three different baryonic 
density values leading to the appearance of two triple points in the phase diagram which have not 
been observed before. Another novel result is that at high pressure and low temperature we observe 
a first order transition which presents a negative slope in the P — T that is reminiscent of the solid- 
liquid transition line observed within the water phase diagram. These unusual patterns occur for 
Gy and B values which lie within the range presently considered in many investigations related to 
the study of magnetars. 

PACS numbers: ll.lO.Wx, 26.60.Kp,21.65.Qr, 25.75.Nq 


I. INTRODUCTION 

Strongly interacting magnetized matter may be produced in non central relativistic heavy ion collisions [1] and may 
also be present in magnetars [2] reaching up about to ^ 10^^ G and ^ 10^^ G in each of these two physical situations. 
As far as heavy ion-collisions are concerned the presence of a strong magnetic field most certainly plays a role despite 
the fact that, in principle, the field intensity should decrease very rapidly lasting for about 1-2 fm/c only [1]. The 
possibility that this short time interval may [3] or may not [4] be affected by conductivity remains under dispute. 

A striking property expected to occur in such extreme conditions is the so called Magnetic Catalysis (MC) phe¬ 
nomenon [5], which implies that the order parameter for the chiral transition represented by the quark-antiquark 
condensate rises as the field becomes more intense (see Ref. [6] for a review). It is then natural to ask how the ex¬ 
pected quantum chromodynamics (QCD) transitions are affected by the presence of intense magnetic fields, a question 
which has been addressed in many recent works (see Refs. [7] for an updated discussion). To summarize the main 
results obtained in these investigations let us start by recalling that at finite temperatures and vanishing chemical 
potential both, model approximations [8-11] and lattice QCD (LQCD) evaluations [12, 13] agree that a cross over, 
which is predicted to occur in the the absence of magnetic fields [14], persists when strong magnetic fields are present. 
However, a source of disagreement between recent LQCD evaluations [13] and model predictions regards the behav¬ 
ior of the pseudocritical temperature (Tpc), at which the cross over takes place, as a function of the magnetic field 
intensity. The LQCD simulations of Ref. [13], performed with 2 + 1 quark flavors and physical pion mass values, 
predict that Tpc should decrease with B while early model evaluations predict an increase (see Ref. [9] and references 
therein). This problem has been recently addressed by different groups [15] which basically agree that the different 
results stem from the fact that most effective models miss back reaction effects (the indirect interaction of gluons and 
B) as well as the QCD asymptotic freedom phenomenon. On the other hand, LQCD results are currently unavailable 
at high densities and low temperatures so that one has to rely in model approximations [9-11] which predict that 
the first order chiral transition takes place at a coexistence chemical potential value which is lower than the one 
observed for the case of unmagnetized matter leading to the phenomenon called Inverse Magnetic Catalysis (IMC) 
[16] (see Ref. [17] for a physically intuitive discussion of the IMC phenomenon). Further investigations performed 
with the Nambu-Jona-Lasinio model (NJL) have been carried out to locate the critical end point (CP) as well as the 
coexistent chemical potentials associated with the first order chiral transition [9, 11]. The results suggest that the 
size of first order transition line increases as the field becomes stronger affecting the position of the (second order) 
CP where the first order transition line terminates. At the same time, the size and location of the coexistence region 
in the presence of a magnetic held appears to oscillate around the B = 0 values [18]. Together, all these effects have 
interesting consequences for quantities which depend on the details of the coexistence region such as the surface ten¬ 
sion and may have consequences for studies related to the properties of magnetized compact stellar objects [19]. The 
motivation for the present investigation stems from an early work by Ebert and collaborators [20] who recognized that 
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the MC phenomenon, associated to the filling of Landau levels, could lead to more exotic phase transition patterns 
as a consequence of the induced magnetic oscillations. To confirm this assumption these authors have considered a 
wide range of the (scalar) coupling values for the two flavor NJL model at vanishing temperatures and, as expected, 
have observed unusual phase structures as a function of the chemical potential such as an infinite number of massless 
chirally symmetric phases, a cascade of massive phases with broken chiral invariance among other features. More 
recently, this seminal study has been extended by a more systematic, and numerically accurate analysis with two and 
three flavors in Refs. [21] and [22] respectively. The results confirm that for certain model parametrizations one is 
able to observe more than one first order phase transition, which is signaled when the thermodynamical potential 
develops two degenerate minima at different values of the coexistence chemical potential ^. It is important to remark 
that, in general, weak first order transitions can be easily missed in a numerical evaluation due to the fact that 
the two degenerate minima appear almost at the same location being separated by a tiny potential barrier so that 
their study requires extra care. Physically, this corresponds to a situation where two different (but almost identical) 
densities coexist at the same chemical potential, temperature and pressure. Here, one of our goals is to extend the 
investigation of these cascades of first order phase transitions, observed in Refs. [20-22], when hot magnetized quark 
matter is subject to the presence of a repulsive vector channel parametrized by Gy. This type of interaction provides 
a saturation mechanism similar to those found in effective nuclear models [24] and is known to be important for an 
accurate description of quark matter at high baryonic densities [25]. As we have verified in a previous work [18], 
the repulsive vector coupling modifies the magnetic effects mainly at lower temperatures and plays an opposite role 
compared to B in the QCD phase diagram. As a matter of fact, the increase of the magnetic field shifts the first order 
transition to lower values of the coexistence chemical potential while a nonzero vector repulsion produces the reverse 
effect. This feature has also been recently verified within the three flavor case [26] in the analysis of compact stellar 
objects. Here, we shall see that the presence of a vector repulsion allows for further interesting possibilities associated 
with the chiral first order transition due to the fact that this term can stabilize intermediary density magnetic phases. 
Being carried out at finite temperatures the present investigation also allows for a more complete description of the 
coexistence region and makes it possible to better understand the physical nature of the associated phase transitions. 
As a first novelty we show that within a cascade of the transitions the one which takes place at the highest pressure 
value seems to be reminiscent of the “solid-liquid” type of transition displayed by the phase diagram of water while 
the others are of the usual “liquid-gas” type (commonly observed within QCD effective theories). Motivated by the 
fact that the Lennard-Jones potential, which describes water, also has a repulsive part we have scanned over the Gy 
values to search for any eventual triple point in the phase diagram of magnetized quark matter since this situation 
cannot be completely ruled out in the scenario considered here. The reason is that the presence of a magnetic held 
induces the free energy to develop multiple minima while the repulsive interaction favors stability so that eventually 
three (instead of the usual two) minima could be globally degenerate leading to the coexistence of three phases. As we 
shall demonstrate the numerical results have confirmed our expectations so that for very particular, but yet realistic, 
parameter values the phase diagram of strongly interacting magnetized quark matter may indeed display triple points. 

The work is organized as follows. In the next section we present the free energy in the presence of a repulsive 
vector channel for the magnetized NJL model within the mean held approximation (MFA) framework. In Sec. HI we 
discuss how a cascade of the usual type of first order phase transitions, with two degenerate minima, takes place at 
finite temperatures. In the same section we extend the analysis to the finite temperature domain in order to draw 
the phase diagrams in the T — /i, T — pB and P — T planes. Section IV is devoted to the study of unusual first order 
phase transitions where the free energy develops more than two degenerate minima allowing for the existence of triple 
points. Our conclusions and final remarks are presented in Sec. V. 


II. NJL MAGNETIZED FREE ENERGY WITH A REPULSIVE VECTOR INTERACTION 

The QCD interaction between quarks can be effectively described by the well-known NJL theory [27] which repro¬ 
duces, at lower energies, the main features of chiral symmetry breaking. In the usual two flavor version the same 
coupling constant, Gs, sets the interaction strength in both the scalar and pseudo-scalar channels. However, in finite 
density investigations the model produces more realistic results if an additional repulsive vector channel, parametrized 
by Gy, is introduced. In this case, the corresponding Lagrangian density [25, 28] can be written as 

- to)V’ + , (2.1) 


^ This fact has also been recently observed to arise within another effective four fermion theory described by the 2 + 1 d Gross-Neveu 
model [23]. 
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where m = rUu — rrid is the bare quark mass. In order to derive the thermodynamical potential within the MFA the 
quadratic interaction terms appearing in the above Lagrangian are linearized by the introduction of the mean fields 
expressed in terms of the scalar condensate, (j) = ('^'0), and the quark number density, p = ('0^'^) 

— (j)^ and 2pip^^^ — , (2.2) 

where quadratic terms in the fluctuations have been neglected while the pseudo-scalar term does not contribute at 
this level. Then, in the case of symmetric quark matter [p = the theory is described by 

^ -M + , (2.3) 

where the effective quark mass, M, and the effective chemical potential, /i, are determined upon applying the cor¬ 
responding minimization conditions, SCl/SM = 0 and SQ/Sp = 0. Integrating over the fermionic fields yields the 
thermodynamical potential 


n = 


{M -mf 

AGs 


AGv 
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d'^p 

(27r)4 


ln[-p2 + m2] . 


(2.4) 


One can then include the effects of a static magnetic field and a thermal bath to this dense quark matter system by 
applying the following replacements [29] to Eq. (2.4): 


po ^ - 'ilA 5 


^ pI + (2n + 1 — s)\qf\B , with s = ±1 , n = 0,1, 2... 
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where ujy = {2ip ^ 1)t:T {u = 0,±1,±2...) represents the Matsubara frequencies for fermions. The Landau levels are 
labelled by n while the absolute values of quark electric charges \qf\ are \qu\ = 2e/3 and \qd\ = e/3 with e = 1/a/ 137 
representing the electron charge Then, following Ref. [30] we can write the thermodynamical potential as 


n = 


(M — m)^ {p — pY 


AGs 


4Gy 


Vt 


vac ~ ^raag 




^med 5 


(2.5) 


where the vacuum contribution reads 


^vac = —2NcNf 


/ 


d^p 

(^ 


+ M2 


( 2 . 6 ) 


and, as usual, can be regularized by a non-covariant sharp cut-ofF, A, yielding 


^vac- g^2 




(A -j- eA) 

M 


— eAA[A^ H- eA^] 


(2.7) 


where eA represents the energy a/A^ H- at the cutoff momentum value A. We remark that in Refs. [31, 32] the 
authors choose a smooth cut off to avoid unphysical oscillations, which appear when the pairing interaction is included 


^ Our results are expressed in Gaussian natural units where 1 MeV^ = 1.44 x 10^^ G. 
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because the sharp cut off limits the allowed momenta. In the present work, no superconducting phase (that would 
require the pairing gap equation to be solved) is used and hence we do not face the problem of unphysical solutions. 
The magnetic contribution of the thermodynamical potential is given by 

^mag = - Y, —^ 2 ^/] - -Xf)lnXf + | , (2.8) 

whre we have used the definition Xf = M‘^/{2\qf\B) and the derivative of the Riemann-Hurwitz zeta function 
= d({z^Xf)/dz\z=-i (see the appendix of Ref. [30] for detailed steps). Finally, the term flrned repre¬ 
sents the in-medium contribution 


^med 


aki\qf\B) j + + , 


f=uk=0 


(2.9) 


where ak = 2 — Sko and Ep = y/^~^^~2E\qJ\B~^^^W^. A similar expression for the magnetized thermodynamical 
potential at Gy = 0 was originally obtained in Ref. [33] where Schwinger’s proper time approach has been used. 
Solving SCl/SM = 0 and SCl/Sjd = 0 we get the following coupled self consistent equations 

M = m-2Gs(l) , (2.10) 

and 

ji = li-2GvP , (2.11) 


Note also that, in principle, one should have two coupled gap equations for the two distinct fiavors: Mu = rriu — 
2Gs{{uu) -h {dd)) and = rrid — 2Gs{{dd) + {uu)) where {uu) and (dd) represent the quark condensates which differ, 
due to the different electric charges. However, in the two fiavor case, the different condensates contribute to Mu and 
Md in a symmetric way and since rriu = rud = m one has Mu = Md = M. The quantities 0 = (j)vac + 4^mag + 
and p appearing in Eqs. (2.10) and (2.11) are given by 


MNcNf 

27r^ 




( 2 . 12 ) 
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where 


(2.14) 


(2.15) 


= 1 + e(g.-M)/^ = i + e(^.+A)/T 

represent, respectively, the Fermi occupation number for quarks and antiquarks. 

At r = 0 the relevant in-medium terms appearing in f2, M and jl can be written as 
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and 


d kf^rnax 

4^med — ^ ^ ^ ^ 
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(2.19) 


where kpik^B) represents the Fermi momentum, kp = — Sf{k, B)‘^ , and Sf(k^B) = ^JM‘^ + 2\qf\kB. The 

upper Landau level (or the nearest integer) is defined by 




,max 


. 2|g/|-B _ 


( 2 . 20 ) 


To obtain numerical results we must now fix the model parameters and, as usual, the cut-off value together the other 
parameters G^, m are chosen to reproduce the phenomenological values [28] for the pion mass ^ 140 MeV), 
the pion decay constant ^ 93MeV) and the quark condensate {{'ip'ipyB ^ 250MeV). Here, we consider the set 
A = 590 MeV, G^A^ = 2.435 and m = 6.0 MeV. 

Fixing Gy poses and additional problem since this quantity should be fixed using the p meson mass which, in 
general, happens to be higher than the maximum energy scale set by A. At present, the vector term coupling Gy 
cannot be determined from experiments and lattice QCD simulations (LQCD) but eventually, the combination of 
neutron star observations and the energy scan of the phase-transition signals at FAIR/NICA may provide us some 
hints on its precise numerical value. Therefore, at the present stage the vector coupling Gy is usually taken to be 
a free parameter whose accepted values he within the range 0.25 Gs- — 0.5 G^ [34-36]. It is worth to point out that 
the explicit use of Gy within the NJL model can be avoided provided that the evaluations be performed beyond the 
large-Vc (or MFA) level. In this case two loop (exchange like) terms bring finite Nc corrections to the pressure such 
as {Gs/Nc)p^ so that the same type of physics can be observed with one less parameter [37]. 


III. CASCADE OF FIRST ORDER CHIRAL TRANSITIONS WITH TWO COEXISTING DENSITIES 

Within the usual first order transition scenario the free energy develops a pair of degenerate global minima defining 
two distinct values for the effective quark masses as well as for the quark number density. Usually, when H = 0, only 
one first order chiral transition takes place so that, for a given temperature, a unique value for the coexistence chemical 
potential exists. However, as already emphasized, the presence of a magnetic field causes an oscillatory behavior which 
may induce more exotic patterns like the appearance of a cascade of first order transitions such as the ones studied 
in Refs. [20-22] where only the case of vanishing temperatures, in the absence of repulsion, has been considered. 
One of the main outcomes of these applications is that, for some phenomenologically acceptable parametrizations, 
several transitions are needed to move from the vacuum phase to the (approximately) chirally symmetric phase. More 
recently, Allen, Pagura and Scoccola [38] have observed a similar situation when considering a a generalized version 
of the two flavor NJL model, which includes a Gy term, but again at vanishing temperatures. In this section we will 
extend the investigation performed at Gy 7 ^ 0 to the finite temperature case to gain extra insight on the effects of B 
and Gy by exploring the phase diagram in planes such as T — //, T — and P — T. This exercise will allow us not 
only to review some of the main aspects related to the appearance of a cascade of first order phase transitions as the 
ones studied in Refs. [20-22, 38] but also to explore more realistic finite temperature situations. 

Let us first recall that, within this model, the order parameter associated with the chiral transition is {{uu) -h {dd))/2 
which, in our case, is directly related to the effective mass as Eq. (2.10) indicates. Therefore, it is instructive to analyze 
the effects of B and Gy over M since this quantity also determines the behavior of the associated EoS. The left panel 
of Eig. 1, which was originally obtained in Ref. [18], displays the effective quark mass, at T = 0, as a function of 
p for Gy = 0.2 Gs- and different values of the magnetic field. The figure indicates that, in the vacuum, the value of 
M tends to increase with B which is in accordance with the magnetic catalysis effect [5]. Also, due to the filling of 
Landau levels, one observes the typical de Haas-van Alphen oscillations which are more pronounced for small values 
of B. Only the segments of the curves where dM/dp < 0 correspond to energetically favored gap equation solutions 
and in the present case (Gy = 0.2 G^) we observe that for eB = 5m^ and eB = 8m^ some of these solutions are 
stable leading to intermediate transitions. Eor example, at eB = 5m^ one sees a first transition at p = 388.55 MeV 
when the mass jumps from 409 MeV to 313 MeV. This is followed by another transition from M = 312.6 MeV to 












6 



FIG. 1: Cascade of first order transitions. Left panel: The effective quark mass, M, at T = 0, as a function of ji for different 
values of eB at Gv = 0.2 Gs-. Right panel: The effective quark mass, M, at T = 0, as a function of for different values of 
Gv at eB — 5m^. In both cases the thick lines represent stable solutions to the gap equation. 


M = 190 MeV at // = 390 MeV. To complete the “cascade” of (three) first order phase transitions one observes a 
final transition from M = 150 MeV to 59MeV at /i = 402.65 MeV. The numerical results illustrated by the figure 
clearly display the effective mass oscillatory behavior showing that at relatively weak magnetic fields we observe many 
oscillations of quark mass values (due to the many Landau levels available). When the field becomes stronger, the 
quantity of oscillations in the quark effective mass is reduced since there are less Landau levels available. The left 
panel also shows that, for a fixed value of Gy, the transition to the chiral phase (lowest mass value) occurs at lower 
chemical potential coexistence values as B increases in accordance with the IMG phenomenon. The right panel shows 
that, for a fixed 5, the transition to the chiral phase occurs at higher chemical potential coexistence values as Gy 
increases (for a complete discussion on these issues the reader is referred to Ref. [18]). For our present purposes it is 
important to note how the insertion of a repulsive vector coupling Gy between quarks brings stability (dMjd[i < 0) 
to the system and occasionally promotes these magnetic intermediary density states to stable ones as the right panel 
shows. 




// [MeV] T [MeV] 

FIG. 2: Quark matter phase diagrams in the T — fi plane (left panel) and in the P — T plane (right panel) for eB = 5.1 and 
Gy = 0.3 Gs. In both figures the continuous thin line represents the case with Gy = 0 for comparison. When Gy = 0.3 Gs 
a cascade of three first order transitions appear at very low temperatures. Each first order transition line terminates at the 
indicated critical point. 
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Figure 2 shows the T — ji and P — T phase diagrams for the case of a field strength^ eB = 5.1 and a vector 
coupling magnitude of Gy = 0.3 which is the value recently suggested by Sugano et al. [36], compared to the 
Gv = 0 case. At temperatures close to zero, one observes a splitting of the Gy = 0 first order transition line into 
three lines occurring at different // values as one could expect from the discussion related to figure 1. As already 
emphasized this exotic scenario may also appear in the Gy = 0 case provided that one uses a different parametrization 
for Gs as discussed in Refs. [20-22]. Here, we have instead adopted a rather canonical value for Gs so that only one 
transition occurs if Gy = 0. This standard choice is well suited to achieve our goal since the role played by the vector 
channel itself can be further highlighted. The left panel of this figure shows that, as expected [39], Gy weakens the 
first order transition and shifts the coexistence chemical potential to high values. Note also that for Gy ^ 0 the third 
transition (thick continuous line) is quickly washed out for temperatures higher than « 4.75 MeV. The right panel 
shows the phase diagram in the physically more intuitive P — T plane. It is interesting to note that the transition 
terminating at CP3, which is associated with the Landau level jump kd = I ^ 2, has a negative AP/AT slope 
(related the Clausius-Clayperon equation) while the ones terminating at GPl and GP2, associated with kd = 0 ^ 0^ 
have a positive slope just like the one usually observed when Gy = 0. Therefore, apart from the usual “liquid-gas” 
type of transition (positive slope) it appears that the combined presence of Gy and B may also induce a transition of 
the “solid-liquid” type observed in the water phase diagram (negative slope) which, as far as we know, has not been 
observed before within QCD motivated models. 

Next, in the left panel of Fig. 3 we show the phase diagram in the T — pB plane which could not be analyzed in 
the previous applications at T = 0 [20-22, 38]. The figure shows that for the chosen Gy and B a values a shrinkage 
of the coexistence regions also takes place. At T = 0 the figure shows a coexistence between the following pairs of 
densities: pB = ^ and pB = 0.8 po (at p = 388.55 MeV); pB = 0.9 po and pB = 1.7 po (at p = 390 MeV); and pB = 2 po 
and pB = 2.55 po (at p = 402.65 MeV) whereas for the Gy = 0 case the dilute phase occurs at vanishing density and 
the dense phase occurs at pB = 3po. As usual the baryonic density is defined as pB = p/3 while po = 0.17fm“^. 
According to Refs [19, 40] one can then expect that the surface tension between the coexistence phases will be lower 
when Gy ^ 0. The right panel of Fig. 3 shows a three dimensional plot displaying the P — T — pB phase diagram 
with the Andrews isotherms that define the equation of state. This figure indicates that the cascade of three first 
order transition observed so far always refer to the coexistence of two phases occurring at distinct pressures. 



FIG. 3: Left Panel: Coexistence phase diagram in the T — ps plane for eB = 5.1 and Gy = 0.3G^. The case Gy = 0 
(thin continuous line) is shown for comparison. Right panel: The phase diagram in the P — T — ps space. In both panels it is 
possible to distinguish the three independent first order phase transitions defined by the coexistence of two distinct densities 
at the same pressure. 

Up to this point we have seen that, due to the LL filling procedure, chiral symmetry restoration may take place via 


The choice of this particular value will become clear in the sequel. 
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successive first order transitions between different magnetized phases. From the free energy perspective this succession 
happens because, depending on the values of the couplings, the presence of a magnetic field may induce the appearance 
of more than two minima so that by varying a control parameter, such as /i, one may observe multiple transitions. 
This scenario can happen either when Gy = 0 and Gs is relatively weak (leading to effective quark masses such as 
M ^ 200 MeV) [20-22] or when Gy 7^ 0 but Gs leads to more standard M values as we have shown. At this point it 
is crucial to note that, with the values of Gy and B considered so far, we have only observed the occurrence of two 
degenerate global minima signaling the usual type of first order phase phase transition for a given value of /i while any 
eventual extra minima remain local. Then, at another chemical potential value a global minimum turns into a local 
one and vice versa leading us to observe a cascade of transitions where only two densities coexist for a given value of 
yU. The results found in this section suggest that perhaps there are certain values of Gy at which three degenerate 
global minima will emerge signaling the coexistence of three different phases at the same T, // and B values as we shall 
discuss next. 


IV. FIRST ORDER CHIRAL TRANSITIONS WITH THREE COEXISTING DENSITIES 

The previous discussion shows that the B and Gy values considered so far produce a cascade of first order transitions 
where, at a given coexistence // value, only a pair of degenerate (global) minima coexist with other (local) minima. 
However, for other parameter values, it may be energetically preferable that one of these local minima becomes a 
global one so that the ground state is triply degenerated. With this motivation let us now scan Gy and eB around 
the values Gy = 0.2 Gs and eB = 5m^. Fig. 4 shows the thermodynamic potential evaluated at eB = 5.1 and 
Gy = 0.2 Gs- at various temperatures. To facilitate the understanding of the unusual phase diagrams which appear 
in the sequel let us discuss the results shown in this figure with some detail. Starting with the T = 0 case one 
observes (left panel) three coincident minima occurring at the same null pressure coexistence chemical potential value 
hTPi = 388.05 MeV where the subscript stands for “triple point 1”. 




FIG. 4: Thermodynamic potential at Gy = 0.2Gs and eB — 5.1 for some selected temperatures. Left panel: at T = 0 
the potential is triply degenerate at the coexistence chemical potential value iitpi — 388.05 MeV. At T — 10 MeV the 
potential is doubly degenerate at two different coexistence chemical potential values, fiia = 387 MeV and fiih = 387.65 MeV. 
At T = 21.4 MeV it becomes triply degenerate again, at iJiTP 2 = 384.82 MeV while at T = 30 MeV it is doubly degenerate at 
l^ii = 381.08 MeV. Right panel: 

At T = 10 MeV this triple coexistence at the same chemical potential no longer survives and one observes the 
dissociation into two separate first order transitions occurring at two distinct coexistence chemical potential values 
hia = 387 MeV and fiib = 387.65 MeV, respectively which signals a “cascade” of two subsequent first order phase 
transitions. At a higher temperature, near Ttp 2 = 21.4 MeV, three degenerate minima are again observed and one 
can ascribe another triple point taking place at iitp 2 = 384.82 MeV. Above this temperature, we observe the usual 
behavior of the chiral with only one first order transition which evolves to the critical point as T increases. The right 
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panel of Fig. 4 shows a final transition line starting at T = 0 and iini = 404.70 MeV and terminating at the CP 
located at Tcp 2 = 4.75 MeV and /icP 2 = 401.60 MeV. 

Notice that in the previous section we have deliberately shown results for the cases Gy =0.2 MeV and eB = b MeV 
as well as for Gy = 0.3 MeV and eB = 5.1 MeV but not for the present choice, Gy = 0.2 MeV and eB = 5.1 MeV, 
which induces the appearence of triple points. As already discussed the quark effective mass is directly related to the 
order parameter so that the discussion can be further clarified by investigating how this quantity varies with /i for 
the T, Gy and B values considered in Fig. 4. With this aim we present Fig. 5 where one may observe the transitions 
discussed in connection with the former figure from a different perspective. The squares locate the first triple point, 
the triangles locate the second triple point while the ordinary type of first order transitions is denoted by the dots. 



FIG. 5: Sequence of plots showing first order chiral transitions for eB = 5.1 and Gy = 0.2 Gs. The squares and triangles 
indicate triply degenerate transitions while the dots indicate doubly degenerate ones. 

From the sequence shown, it becomes evident how the triple point TPl, related to three different mass values 
(marked by the squares), determines the same chemical potential fippi = 388.05MeV when T = 0. Note that the 
the three mass values, satisfying the gap equation, differ by an approximately equal amount (close to 100MeV). 
Still at T = 0 but at a higher chemical potential value one observes the occurrence of another (ordinary) transition 
signaled by two different mass values (marked by the dots) at jiin = 404.70 MeV. Then, rising the temperature to 
T = 10 MeV one sees that the triple coexistence, observed at T = 0, decouples into two ordinary first order transitions 
occurring at two different, but yet very similar, chemical potentials given by /ija = 387 MeV and fipj = 387.65 MeV. 
At this temperature the first order line starting at /i/// and T = 0 has already vanished. Note also that the low mass 
value occurring at fiia and the high mass value occurring are almost identical (the difference amounts to about 
5 MeV). Near T = 21.4MeV, these distinct first order lines start to approach each other and one observes another 
triple point (TP2) at iitp 2 marked by the triangles. In this case the high and the intermediate mass values become 
almost identical but differ substantially from the low value (by almost 200MeV). Above this temperature, as the 
T = 30 MeV panel suggests, only the ordinary scenario takes place until the first order transition line ends at a critical 
point, CPI, located at Tcpi = 57.65 MeV and ficpi = 361.90 MeV. 
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We are now in position to map all these transitions into phase diagrams such as the ones shown in Fig. 6 in order 
to display the phase boundaries in the T — /i and P — T planes. 




FIG. 6: Phase diagrams in the T — /x plane (Left Panel) and in the P — T plane (Right Panel) fovB = 5.1 and Gv = 0.2 Gs- 
At T = 0, both diagrams show a triple point (TPl) splits into two common dual-phase coexistence lines as the zoomed areas 
show. At Ttp 2 = 21.4 MeV, the dual coexistence lines converge again and merge at another triple point (TP2). Beyond this 
temperature one observes the ordinary first order line which ends in a critical point (CPI). The figure also shows that the first 
order phase transition line which terminates at CP2 has a negative AP/AT slope. 

Analyzing these figures we verify that the usual phase diagrams are drastically modified from the ordinary case at 
Gv = 0 since we now have a first-order coexistence line originating from a triple point at (T = 0, jitpi = 388.05 MeV) 
and then splitting into two new branches of ordinary first order transitions, as shown in the zoomed region, to finally 
merge again at a second triple point {Ttp2 = 21.4 MeV, /iTP2 = 384.82 MeV). Then, for higher temperatures the 
transition always follows the usual pattern until the line ends a critical point (CPI). The figure also shows an 
ordinary first order phase transition line which starting at (T = 0,/i/// = 404.70 MeV) and terminating at CP2 
{Tcp2 = 4.75 MeV, jiicP2 = 401.60 MeV). The P — T allows to further understand these two different first order lines 
by observing the dP/dT slope in each case. The diagram shows that the line terminating at CPI has a positive slope 
just like in the well known “liquid-gas” type of transition while the line terminating at CP2 presents a negative slope 
which is reminiscent of the “solid-liquid” transition occuring in the water phase diagram which does not appear to 
have been reported previously (at least within the context of strongly interacting systems). The left panel of Fig. 7 
displays the coexistence phase diagram in the T — pB plane while the right panel of the same figure shows a three 
dimensional plot displaying the P — T — pB phase diagram with the Andrews isotherms. 

It is well known that at zero temperature this type of model, as well as the quark meson model, typically predict 
phase coexistence between the vacuum and the dense quark phase, similar to the “liquid-gas” transition. However, 
we now observe that when a vector repulsion and a magnetic field are present it is possible that three phases coexist 
for particular parameter values. As Fig. 7 shows, for eB = 5.1 ,m^ and Gy = 0.2Gs^ one observes the coexistence 
between phases with densities Pb/P o = 0, 0.85 , and 1.75 which is a rather interesting result since the vanishing density 
represents the vacuum, the intermediate density value is very close to that of ordinary nuclear matter while the third 
is close to the values this model predicts for the dense quark phase. Eventually, by scanning over parameter values one 
may force the intermediate density to take place at = 1 so that this could be identified as nuclear matter. As the 

temperature increases, the three-phase coexistence vanishes and one observes a decoupling into two separate branches 
la and Ib^ which are related to two first order transitions taking place at distinct pressure values. This bifurcation 
in two coexistence branches creates a stability island between the regions la and Ib. Again, at Ttp 2 = 21.4 MeV 
these two distinct first order transitions will share the same pressure at a triple point TP2 as the figures show. The 
remaining first order transition corresponding to the coexistence region II ends at an ordinary critical point GPl. 
Furthermore, as already emphasized, another coexistence region appears at a higher density range, labelled as III in 
the left panel of Fig. 7. At T = 0 the coxistence densities for this region are Pb/P o = 2.0 and Pb/P o = 2.6. Then, 
when the temperature reaches the critical value Tcp 2 = 4.75 MeV this coexistence region, which corresponds to a 
“solid-liquid” type of transition, terminates at GP2. 

At this point it is important to recall that the occurrence of multiple phases at the same pressure is only possible if 
thermodynamical potential develops multiple global minima. This condition can be attained to a particular choice of 
B and Gy values which give rise to multiple densities coexisting in the same transition chemical potential p. Let us 
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FIG. 7: Phase diagram in the T — ps plane (Left Panel) and EOS isotherms in the P — T — ps space (Right Panel) for the 
B = 5.1 and Gv = 0.2 Gs- case. Filled squares and triangles indicate the triple points and filled circles indicate the critical 
end points in accordance with the previous figures. 


now call Gy the critical vector coupling value at which three degenerate global minima appear for a given B value. 
Then, when mapping the parameter space of (5, Gy, p) values which result in a three-phase coexistence at T = 0 we 
encounter the oscillatory behavior shown in Fig. 8. 



FIG. 8: Magnetic field dependence of the critical values, Gy, related to the transitions (ku, ka) (0,0) and (/cn, kd) (0,1) 
at T = 0. The dot labels the values used in our analyzes and the open circles mark the excluded Gv = 0 case. 

Each critical vector coupling Gy can be associated with the emergence of a new intermediary (stable) value of the 
order parameter. Depending on the field strength, more than one intermediary Landau level configuration {ku^ kd) 
can appear and each one has its own Gy value. So, it must exist an infinite set of critical values Gy associated with 
the infinite number of intermediary states {ku, kd) in the limit of B ^ 0 and Gy oo. The figure only shows the 
first two possibilities characterized by the transitions (/c^, kd) (O^O) and (/c^, kd) As one can observe 

the former transition covers the region of the parameter values which are usually considered in the recent literature 
(Gy = 0 — 0.5 Gs and eB = 2 —8m^). Note that for high fields eB ^ 9m^ the emergence of triple points is suppressed 
since only the LLL is occupied. We have tested coupling values above and below the critical line observing that when 
Gy < Gy only the traditional case with one first order transition line emerges. Then, when Gy = Gy the multiple 
coexistence becomes possible and a new first order transition lines starts from the original one at T = 0 as shown 
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in Figs. 6 and 7. Finally, when Gy > Gy the chiral symmetry restoration takes place via more than one first order 
transition line in a cascade mode as discussed in Section III. We close this section by remarking that we have also 
observed the coexistence of four phases but this is resctrited to the T = 0 case only. 


V. CONCLUSIONS 

We have investigated how the presence of a magnetic field and a repulsive vector interaction may influence the 
phase diagram of strongly interacting matter generating unusual transition patterns associated to first order chiral 
transitions. In the first part of the present study we have considered parameter values which produce a cascade of 
first order transitions similar to the ones analyzed in Refs. [20-22] at vanishing temperatures. We have taken a step 
forward by incorporating the vector interaction as well as by pushing the evaluations to the finite temperature domain 
in order to better understand the physical nature of such transitions. Mapping the transition into the P — T plane 
has allowed us to observe that the transition which takes place at high pressure has a negative dP/dT slope which is 
reminiscent of the “solid-liquid” transition observed in the water phase diagram while the remaining (low pressure) 
transitions have a positive slope and therefore resemble the “liquid-gas” transition which is usually observed within 
effective quark models. Having in mind that, due to the inherent oscillations caused by the filling of Landau levels, 
the magnetic field induces the free energy to develop multiple minima while the repulsive vector interaction favors 
stability we have scanned over B and Gy values to check for the possibility of observing three degenerate global 
minima which would then lead to the existence of a triple point. Our expectation was confirmed by the numerical 
investigation and we were able to find, for a given value of H, a certain value of Gy so that three (instead of the usual 
two) phases coexist which, as far as we know, has not been observed before. In contrast to the “solid-liquid” type of 
transition, which may be find when canonical parametrizations are used, the existence of triple points is only possible 
when very specific values of Gy are chosen. Finally, we point out remark that the coexistence of multiple phases has 
also been recently observed by Manso and Ramos [41] within the 2 1 d Gross-Neveu model despite the fact that a 

repulsive term has not been considered. Instead the authors have considered a tilted magnetic field observing that the 
parallel component tends to stabilize the free energy (just as Gy in our case) so that three degenerate minima were 
also observed. Together with our observations this result leads us to conclude that the phase diagram of magnetized 
matter may display the coexistence of multiple phases if the dynamics contains terms which bring stability to the free 
energy. 
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